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Abstract: We consider the moduli space of 1/2 BPS configurations of type IIB SUGRA 
found by Lin, Lunin and Maldacena (hep-th/0409174), and quantize it directly from the 
supergravity action, around any point in the moduli space. This quantization is done 
using the Crnkovic-Witten-Zuckerman covariant method. We make some remarks on the 
applicability and validity of this general on-shell quantization method. We then obtain 
an expression for the symplectic form on the moduli space of LLM configurations, and 
show that it exactly coincides with the one expected from the dual fermion picture. This 
equivalence is shown for any shape and topology of the droplets and for any number of 
droplets. This work therefore generalizes the previous work (hep-th/0505079) and resolves 
the puzzle encountered there. 
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1. Introduction 



One of the greatest achievements of string theory is the fact that it is able to explain 
gravity, or supergravity, in terms of a microscopic set of degrees of freedom — oscillation 
modes of strings with various boundary conditions. 

Undoubtedly, if string theory is to be a quantum theory of gravity, then such a de- 
scription is essential. One particular application of having a microscopic set of degrees of 
freedom amenable to quantization is the ability to count the number of microstates cor- 
responding to a specific geometry, thus determining its microscopic entropy. One of the 
most fascinating applications of this is for geometries describing black holes. 

In general a coherent state of string oscillations generates a supergravity background. 
However, there are very few such backgrounds for which we know the corresponding co- 
herent state of string oscillations. In most cases, given some supergravity background, it 
will be very hard to find the exact coherent state making it up. In some specific cases, 
when the spacetime is asymptotically AdS, one could use AdS/CFT Q to obtain a different 
description of the supergravity configuration — in terms of microscopic degrees of freedom 
living in a conformal field theory, or in a deformation thereof. Such a description requires 
a good knowledge of the boundary theory and of the bulk/boundary dictionary, which is 
not necessarily available. 

It is therefore clear that it would be very useful if one could quantize geometries directly 
from the supergravity picture, without recourse to a different microscopic description of the 
system. In this paper, continuing the previous work done in collaboration with L. Grant, 
J. Marsano and K. Papadodimas in we develop a method to do precisely that, which 
we propose to call on-shell quantization. 

As the name suggests, on-shell quantization is not a method to quantize all possible 
fluctuation modes of supergravity, but is used only to quantize a subspace of fluctuations 
all describing on-shell SUGRA configurations within a given moduli space. In this sense it 
is a special case of mini-superspace quantization The method is only applicable when 
on the moduli space of solutions all fluxes are kept flxed. 

When applicable, this method is used to quantize only the moduli space of solutions, 
freezing all other fluctuations of the flelds. Thus one might wonder how meaningful the 
results of such a quantization are in the general framework. In general indeed one should 
expect corrections to the results of the on-shell quantization, and one can try to estimate 
their size. It is only in particular cases where the moduli space of solutions is also protected 
by some symmetry, that this sector could completely decouple from the rest of the Hilbert 
space, and the energy spectrum computed by the on-shell quantization method would 
remain uncorrected. 

In this paper we apply the on-shell quantization method to a certain moduli space of 
Type IIB SUGRA, which both has flxed fluxes and is also protected by supersymmetry, so 
that the results of the on-shell quantization remain uncorrected. This is the set of 1/2 BPS 
solutions found by Lin, Lunin and Maldacena in [Q]. The LLM solutions are in one-to-one 
correspondence with various 2-colorings of a 2-plane, usually referred to as 'droplets'. The 
deformations of the droplet shapes such that the areas of different white and black regions 
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remain fixed are exactly the deformations which leave the fluxes fixed. Thus what we would 
like to do is to quantize the moduli space of such fluctuations around a general droplet 
conflguration^ . 

The LLM solutions have the feature that if the droplets are of finite size, then the 
spacetime is asymptotically AdS^ x S^. As we remarked before, this asymptotic behavior 
enables one to use AdS / CFT and get an equivalent description in terms of the = 4 SYM 
field theory. The specific sector of 1/2 BPS solutions considered by LLM turns out to be a 
subsector of SYM, admitting description in terms of free fermions in a harmonic oscillator 
potential 0, |8|. In fact the droplet describing the supergravity solution turned out to be 
the same droplet in the phase space describing the free fermions. Thus this is one of the 
special cases where string theory provides us with a dual description of the system where 
microscopic degrees of freedom are manifest. It is the Hilbert space of these free fermions 
in a harmonic oscillator, which we would like to derive entirely from the SUGRA picture 
using the on-shell quantization. 

In a previous paper with L. Grant, J. Marsano and K. Papadodimas we gave general 
expressions for the symplectic form on the moduli space of LLM solutions, and applied them 
to two specific backgrounds: AdS^ x and the plane wave. In the case of AdS^ x we 
completely reproduced the fermion symplectic form from the SUGRA symplectic currents. 
In the plane wave case we found the symplectic form of the same functional expression 
as we expect from the fermion analysis, however the numerical coefficient was a factor of 
2 smaller than expected. In that paper we gave a few speculations on the source of that 
discrepancy. 

In this paper we resolve the puzzle. We make the observation that one must work in 
a gauge where the variations of the SUGRA gauge fields are everywhere regular. Different 
gauge choices result in extra boundary terms in the symplectic current, and thus working 
in a singular gauge, as we did in might result in missing terms in the symplectic form. 
In this paper we succeed in evaluating the SUGRA symplectic form around any of the LLM 
backgrounds, finding the expected fermion symplectic form with exact numerical matching. 

This paper is organized as follows. In section 2 we first discuss the general issue of 
quantizing geometries from the supergravity action. We explain the motivation and main 
ideas of the 'on-shell quantization' method, and show how to apply it to pure gravity as 
well as to general Lagrangian theories. A central role in this discussion is played by the 
CWZ symplectic currents, introduced 20 years ago by Crnkovic and Witten and by 
Zuckerman We then turn to discuss when this method is applicable, and how one 

may expect it to be corrected by quantization of the full theory. 

Then in section 3 we apply the general idea and expressions to the case of the " Bubbling 
AdS" geometries found by LLM Q. We first analyze the dual fermion system and find an 
expression for the symplectic form in the large limit. This is the expression we aim to 

^Note that we strive to have a quantization which relies only on SUGRA and not on some microscopic 
string or brane picture, such as was done for instance in or, in a different context, in It should be 
noted that although some remarks on the problem of performing a SUGRA quantization were made in 
they seem to apply only in a partial case when there is a preferred group action on the moduli space of 
solutions. In this paper we follow a different, more direct, path. 
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independently derive from the supergravity analysis. Then we indeed turn to the moduli 
space of supergravity solutions, and show that using the CWZ currents we are able to 
reproduce this form, first in the specific case of the plane wave geometry, and then in the 
most general droplet case. We end in section 4 with conclusions and directions for future 
research. Some technical details are deferred to the appendix. 

Results of this paper were reported in a talk presented by the second author at Strings 
2005, Toronto. 



2. Geometry quantization from supergravity 
2.1 Motivation 

Given a continuous family of supergravity solutions, it is natural to ask how the moduli 
space of this family is going to be quantized. A well-understood special case of such 
quantization is the Dirac quantization condition imposed on the fluxes of gauge fields 
present in the theory: typically, such fluxes T on closed cycles have to be quantized in 
units of an elementary flux: T = J-qu, n G Z. 

However, sometimes the moduli space will contain deformations corresponding to all 
fluxes kept fixed. In this paper we are mostly concerned with how to quantize the moduli 
space in such a situation, which is in a sense complementary to the flux quantization. 

As a characteristic example, consider the family of D1-D5 solutions with angular mo- 



mentum found in |11]. The moduli space of solutions is parametrized by 4-dimensional 
closed curves; there are no nontrivial fluxes. This family plays an important role in 
Mathur's program of describing black hole microstates by regular geometries, providing 
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microstates for the 2-charge D1-D5 black hole |12]' 

Another interesting recent example is the LLM family of Type IIB solutions with 
AdS^ X asymptotics |^^. The moduli space is parametrized by planar droplets of 
various shapes. There is an infinite-dimensional family of deformations within the moduli 
space corresponding to keeping the fluxes fixed: these are the deformations keeping fixed 
areas of all black and white regions. We will discuss these solutions and their moduli space 
quantization in detail in Section |3|. 

For both of the above families, there is a dual description — free fermions in the LLM 
case, chiral fundamental string excitations in the D1-D5 case — which can be used to 
quantize the moduli space. These dual descriptions can be derived by using AdS/CFT or 
various other indirect methods (giant graviton picture for LLM, string dualities for BI- 
DS). However, as we will explain below, there exists a method to derive such quantization 
results directly from supergravity. In some cases this might be the only way to quantize 
SUGRA systems, since a dual microscopic description is not always known. 



^In jlsj another family of SUGRA solutions is constructed, which is supposed to represent additional 
microstate geometries. 

'^LLM Q also found similar M-theory solutions with AdSi x and AdSr x asymptotics, and other 
families of solutions have been generated subsequently by applying the LLM method to different theories 
in various dimensions jlj] , , [|l^ . These families became collectively known as "Bubbling AdS". 
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Counting 3-charge D1-D5-P black hole microstate geometries in the context of Mathur's 
program could provide an example of a situation when our direct method may lead to 
progress which will otherwise be difficult to achieve. So far only some specific families 
of regular 3-charge geometries were constructed [17|. It is conjectured that the general 



case can be understood in terms of supertubes [^], however a world- volume description 
of the 3-charge supertube configurations, which could be used to quantize them, is not yet 
known^. In that case our direct method may be the only way to quantize and count these 
geometries, once they are fully described. 

2.2 General idea 

The general idea of the method is simple. Every supergravity theory, as any Lagrangian 
theory, comes equipped with a symplectic form 17, which is defined on the full phase space 
of this theory. The given moduli space of solutions Ai, which we would like to quantize, 
forms a subspace of the full phase space. All we have to do is to restrict Q to A4, which 
will give us the symplectic form on M: 

uj = n\M- (2.1) 

Once u; is found, it can be quantized in the usual way. We will call this method on-shell 
quantization^ . 

This philosophy is completely general and applies to any Lagrangian theory. For 
example, let us quantize the chiral sector of a 2-dimensional free boson using this method. 
The action is 

S = - dtdx(^'^ -(t)'^) , (2.2) 
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and the symplectic form is given by the familiar expression 



n= dx6TT{t,x) A5(f){t,x), TT = (j). (2.3) 

J t=const 

This is a 2-form on the full phase space of the theory (which can be thought of as the 
space of solutions of equations of motion). The 1-forms 5(1) and 5tt (where 6 is the exterior 
derivative operation) can be thought to represent fluctuations around the solutions. We can 
either use the wedge product A as in (2.3) or think of 1-forms as anticommuting quantities. 



Now we should pick the moduli space of solutions to be quantized, which we choose to 
be the space of left-moving solutions, parametrized by an arbitrary function: 

M = {(t){t, x) = f{t + x), f arbitrary}. (2.4) 

It is a simple matter to restrict il. to M. We get: 

u; = n\M = [ du5f'{u) A6f{u). (2.5) 



^There has also been some other work trying to identify CFT states with particular 3-charge D1-D5-P 
geometries [ p^ . 

^ We called it 'minisuperspace quantization' in [Q, however this term is used in a slightly wider sense 
in the canonical gravity literature, implying a reduction of degrees of freedom (typically by imposing a 
symmetry), but not necessarily explicit knowledge of all solutions. 
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Now from to we can infer the Poisson bracket: 



{f{ui)J'{u2)} = 5{ui-U2). (2.6) 

This Poisson bracket can then be promoted to a quantum commutator using the Dirac 
prescription: 

[f{ui)J'{u2)\=ih5{u^-U2). (2.7) 
The Fock space representation of this commutator is then constructed in the usual way: 

f{u)= / Ji e-'P^ap + h.c, (2.8) 



27r^ 

[ap,al,] = 2Trh5{p-p'). (2.9) 

This way we recover the standard quantization of the chiral boson sector, which is equivalent 
to quantizing the full theory and then putting all right-movers in their vacuum state. Notice 
that the only piece of off-shell information we are using is the full symplectic form (^.31). 
All the remaining computations are done on shell, that is in A^. 

2.3 Symplectic form of gravity 

We are interested in theories which contain gravity as a subsector. Thus a necessary 
prerequisite is the symplectic form of pure gravity, which we will now discuss. We will 
assume that the solutions we have to quantize are regular, and that an initial value surface, 
S, can be chosen. We will also assume that there are no horizons present^. All these 
assumptions will be satisfied in the examples to be considered below. 

The most direct route to the symplectic form of gravity lies via the canonical formalism. 



Gravity can be put in the canonical form using the ADM splitting |20| of the metric: 



ds"^ = -N^dt"^ + hij{dx' + N'dt){dx^ + N^dt). (2.10) 

The phase space is parametrized by the metric on S, hij, and by the corresponding canon- 
ically conjugate momenta II*-'. These momenta are related to the extrinsic curvature K^^ 
of the Cauchy surface S: 

U'^ = Vh{K'^ - K\h^i) , (2.11) 
Kij = ^ {hij - D,Nj - DjN,) . (2.12) 

The canonical variables satisfy a set of nonlinear constraints, which have to be realized as 
operator relations in quantum theory. This is a highly nontrivial task when quantizing the 
full theory. Fortunately, in our case this won't be a problem: we are quantizing families of 
classical solutions, and all constraints will be automatically satisfied. 



®In presence of horizons one would have to consider a Cauchy surface in the extended black hole space- 
time, which may or may not exist, and could contain additional causally disconnected asymptotic regions. 
The nature of degrees of freedom which are being quantized in this case is not clear to us. Alternatively, 
one could try to treat horizons as boundaries. Additional analysis is required to decide which approach is 
correct. For now we prefer to exclude this case from discussion. 
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In this formalism, the symplectic form is given by the natural expression: 

n = [ d^-^x5W^{t,x)A5hijit,x). (2.13) 

To restrict this symplectic form to the moduli space of solutions Ai, one first has to put all 
metrics of the family in the ADM form, and evaluate the canonical momenta using ( 2.11| ). 



This is doable in principle, but can be rather cumbersome in practice. Fortunately, an 
equivalent method exists, which avoids using the ADM split. 

2.4 Covariant approach 

In the equivalent covariant approach, which is computationally much simpler than the 
direct method outlined in the previous subsection, one expresses the symplectic form as an 
integral of a symplectic current over the Cauchy surface S: 

= y d^if . (2.14) 

The symplectic current of gravity was found by Crnkovic and Witten and is given by 
the following covariant expression containing variations of the metric and the Christoffel 
symbols'^: 

f = -<5rL A 5(^r") + 5r:;„ a 5(^5'-) . (2.15) 

This current has a number of remarkable properties. First of all, it is conserved: 5; = 0, 
which makes invariant under variations of S (local variations in general, as well as 
global variations if the metric perturbations have fast enough decay at infinity). Second, 
changes by a total derivative if the metric perturbation undergoes a linear gauge trans- 
formation: 

Sgmn Sgmn + ^ (m^n) ■ (2-16) 

This property renders Q defined by ( |2.14| ) gauge invariant (assuming that has finite sup- 
port or decays sufficiently fast at infinity) . It should be noted that both properties are true 
only if the metric is varied inside a moduli space of solutions, i.e. the background satisfies 
Einstein's equations, while 6gmn solves the linearized equations. In general, evaluating the 
symplectic form away from the solution space would be meaningless. 

2.5 Generalization to any Lagrangian theory 

The above strategy admits a natural generalization to any Lagrangian theory @,[10| 



(see also [^). For our purposes it will be enough to assume that the Lagrangian L = 
L^cf)"^ , di(f)^) (where the index A numbers the fields) does not contain second- and higher- 
order derivatives. Under these conditions, the Crnkovic- Witten- Zuckerman symplectic 
current is defined by 



^Notice that Q defines the symplectic form as / dEiy^ —gj\ so that our symplectic currents differ from 
by a factor of y/—g. 
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The symplectic form is defined by ( 2.141) , as before. Both properties — S-independence and 



gauge invariance — are still true, provided that the equations of motion are satisfied. The 
gravitational current (2.15) is obtained from the general formula ( 2.17D , provided that one 



drops the total second derivative term from the Einstein-Hilbert action, which is equivalent 
to adding the Gibbons-Hawking boundary term |22p. 



Another important invariance property of the symplectic current ( 2.171 ) is that its 
definition is independent of the choice of fundamental fields (f)^. For example, in the case 
of pure gravity we may choose {(j)^} = {gmn} or {(f)^} = {5"*"}, the resulting symplectic 
currents being identical. The precise general statement is that the symplectic current 
(2.17) is invariant under point transformations (i.e. transformations which do not involve 



derivatives of the fields) (jj ^ cj)' = ^{(j)). For a proof it is convenient to rewrite the 
symplectic current as 

dL 

fcwz = Sl\ l' = QQ^Scp\ (2.18) 

and show the invariance of Under a point transformation, 6(p^ transforms contravari- 
antly in the A index, i.e. it is multiplied by the Jacobian: 

84>^^J^bH''. J^b = -^. (2.19) 

The derivatives di4>^ also transform contravariantly. Thus dL/ddi(j)^ will transform co- 
variantly, and the product P is indeed invariant. 

This invariance can be used to demonstrate the equivalence between the pure gravity 
canonical symplectic form ( |2.13| ) and the covariant expression ( p.l5| ). We just have to show 
that the gravitational symplectic form evaluated using ( |2.17| ) reduces to ( p. 131 ) if the ADM 
variables N,Ni,hij are used as a set of fields {(t>^} parametrizing the metric. In a sense, 
this is obvious, because the II*-' given by ( 2.1l| ) were in fact defined by ADM ||2^ exactly 



as the variation w.r.t. hij of the gravitational Lagrangian, which written in these variables 
takes the form 

Ladm = NVh{KijK'^ -K^ + i?^^)) . (2.20) 



In addition, time derivatives of and Aj are absent from ( |2.20 ). Thus we see that the 



integrand of ( 2. 13] ) is nothing but J*, which is the only needed component of the symplectic 
current, since T, = {t = const} in this parametrization^ . 

2.6 Restriction of fixed fluxes 

As we stressed in the beginning, the on-shell quantization method is designed to be used 
in the situations when all fluxes are fixed. Actually more is true - the method, at least in 
the presented form, can be used only in such situations. In other words, the regimes of 
applicability of this method and of flux quantization are mutually excluding. 



*It is convenient to use the explicit form of the Gibbons-Hawking Lagrangian Lgh ~ \/—gg^'' (rSTim ^ 
T\fJ^Yrn) (see e.g. |^^) when performing the computation |^]. 

^Strictly speaking, Ladm 7^ icH, because the J?*-^' term in ( ^.20| ) still contains second-order derivatives 
in spatial directions. However, this difference does not affect the time component of the symplectic current. 
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Figure 1: An example of a 3-dimensional moduli space M foliated by 2-dimensional symplectic 
sheets corresponding to a fixed value of a gauge field flux J- . On-shcU quantization can determine 
symplectic structure on the sheets, but cannot be used to quantize the flux. 



To see this in a simple example, suppose that we use the on-shell quantization method 
to quantize the Abelian gauge field on a 4-manifold of nontrivial topology. From the action 

S = -\j d^x^FrnuF'^'' (2.21) 

we derive the symplectic current using ( p. 17 ): 



f = -6iV^F'n^SA^. (2.22) 

We see that the gauge potential variation, SAm, appears in this expression. In order to 
compute the symplectic form, we have to be able to choose a gauge in which 6Am{x) is 
regular everywhere on the Cauchy surface S. This would be impossible unless the field 
strength variation 5F = d{5A) has zero flux on all closed 2-cycles. It is easy to see that 
this restriction is completely general and always appears when there are gauge fields in the 
theory. 

The restriction of fixed fiuxes can be thought of as defining symplectic sheets inside 



the moduli space, which thus becomes a Poisson manifold |24|, with fluxes having trivial 
Poisson brackets with everything else. Quantization inside each of these sheets is governed 
by a corresponding symplectic form. Changing fluxes corresponds to continuously changing 
the sheet. It is only when flux quantization is taken into account, that this continuous 
variation becomes discrete (see Fig. |l|). 

2.7 The role of dynamics 

To get a nontrivial moduli space quantization, the symplectic form computed by the on- 
shell quantization method should be non-degenerate when restricted on M. This implies 
that the canonical momenta should be nontrivial functions on A4, much like in the chiral 



boson example considered in Section 2.2, where we had vr = f'{u). In the case of gravity this 
means that the class of solutions cannot be static. Indeed, for static solutions the extrinsic 
curvatures ( |2.12| ) are zero. Thus the canonical momenta ( 2.11 ) will be identically zero on 



M, and the restricted symplectic form will vanish. This means that for static solutions, the 
parameters characterizing the moduli space do not acquire nontrivial commutators among 
themselves upon quantization. 
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The symplectic form will be usually non-degenerate if some dynamics is present. For 
example, the "Bubbling AdS" solutions considered in Section ^ below have nonzero angular 
momentum: they are stationary but not static, and this will be enough to make the 
symplectic form nondegenerate. The reason is that A'j 7^ in ( p. 12 ) and the extrinsic 



curvatures no longer vanish identically. In [g], the appearance of a nontrivial symplectic 
form on the moduli space was conjectured to be related to super symmetry. However, it 
is easy to see that supersymmetry is neither necessary nor sufficient for that: the moduli 
spaces of cylindrical |2^] or plane [^] gravitational waves have a nontrivial symplectic 
structure in pure gravity; on the other hand, the symplectic form will vanish on any 
moduli space of static supersymmetric solutions. Supersymmetry may be instrumental 
though to control corrections to the on-shell quantization results (see Sections 2.10| ) or to 



obtain simplified symplectic current expressions (see remarks at the end of Section |3.7.3 ). 



The above remarks are useful in order to understand the relation of our method to 



the perhaps more familiar 'moduli space approximation' a la Manton |27, 25, |30| , 
where one usually starts with a moduli space describing multi-centered static black hole 
or soliton solutions. This moduli space is typically finite-dimensional, parametrized by 
the coordinates of the centers and possibly some extra parameters. For example, in the 



multi-centered extremal black hole case [30| the moduli space is 3n-dimensional, where n 
is the number of black holes. Then one considers slow scattering of these black holes (or 
solitons) and finds that it can be described by the geodesies in a certain (calculable) metric 

ds^ = g^j{q)dq'dqK (2.23) 

Then one can consider quantization of this system of slowly moving solitons. But notice 
that in the process of introducing slow motion we enlarged the phase space by adjoining 
canonical momenta (corresponding to velocities). So, in the black hole case, the phase space 
becomes 6n-dimensional, twice the dimensionality of the original moduli space. Also, in 
order to construct these slowly moving black hole solutions (and derive metric ( |2.23| )), 
one should go off-shell, i.e. out of the original moduli space. If one were to recover the 
metric ( 2.23| ) using our method (which may or may not be possible in practice), one would 



have to compute the symplectic form on the extended 6n-dimensional phase space. The 
symplectic form would be degenerate if restricted to the original 3n-dimensional moduli 
space corresponding to the static configurations, in agreement with the above discussion. 
This means in particular that black hole center coordinates are not quantized by themselves. 

2.8 On-shell quantization vs effective action 

In essence, the on-shell quantization method does nothing but quantizing first-order fiuc- 
tuations along the moduli space M (except that it provides a particularly efficient way to 
do this). This is particularly clear from the symplectic current expression ( p. 17 ), which 



involves only first-order perturbations of the fields and, moreover, depends only on the 
quadratic part of the Lagrangian as expanded around a given background. To see this 
explicitly, let us parametrize the fields (j)^ in ( |2.17| ) by their deviations "0"^ from a given 
background 4>q & M.: 

</> = 00 + ^ , (2.24) 
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at the same time expanding the Lagrangian of the theory in terms of ip: 



(2.25) 



1=2 



where L^*) is degree i in i/^. The part hnear in ip is absent, since we assume that (/)o is a 
solution. The coefficients of L^*) may and will generically depend on (pQ, but the precise 
form of this dependence is irrelevant for our argument. 



Now, as we showed in Section 2.5, the symplectic current is independent of the field 



choice. In particular, we can use to compute it. Thus we will have: 

ddiip 



'^cwz = E^j^^SV- (2.26) 

i=2 



bmce is de gree i in ■0, dL^^^ / ddiip will be degree {i — 1), and 6{dL^^^ /dditp) will be 



degree {i — 2) in %p (and linear in 50). So we see that once we evaluate (|2.26) on shell 



by putting ip = 0, all the terms with i > 3 will vanish, and thus the symplectic current 
depends only on the quadratic part of the fluctuation Lagrangian. 

This remark makes it clear that any result obtained by the on-shell quantization 
method can be in principle obtained by a more conventional method based on computing 
the quadratic action, quantizing it, and specializing to a particular subclass of perturbations 
corresponding to deformations along M. In practice, however, there may be formidable 
difficulties in following the conventional path. The main difficulty is that the full quadratic 
action will typically have a complicated form, except around very special backgrounds. 
One could think that perhaps one could consider a truncation of the quadratic action to 
the fluctuation modes along M, and that would have a chance to be simple. The problem 
here is that such a truncation does not even have to exist, since these modes can couple to 
other modes already on the quadratic level. 

These difficulties are best demonstrated on a concrete example of the LLM family of 
solutions Q. On-shell quantization of this family is considered in Section 3 below, and was 
also the subject of Q. In this case the quadratic effective action is available around only 
one representative of the family, namely AdS^ x |31, ^ (and can also be found around 
the plane wave, using the results of |^3|). In any other case it is highly unlikely that a 
tractable quadratic action can be found. Moreover, even around AdS^ x the quadratic 
action couples the moduli space modes characterized by the BPS condition J = A to 
the opposite angular momentum modes J = —A, as a simple consequence of angular 
momentum conservation. This shows that on-shell quantization is the only practical way 
to quantize in the general case. However, around AdS^ x both on-shell quantization and 
the quadratic effective action method can be applied, and the results agree as they should 
i- 

2.9 Hilbert space. Semiclassical states. Hamiltonian 

Using the on-shell quantization method, we can compute the symplectic form on the moduli 
space of solutions. This symplectic form encodes the Poisson brackets between the functions 
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defining the geometry (e.g. the shape of the droplets in the LLM case). Quantizing the 
brackets, we find commutation relations between these functions. We can then find a 
Hilbert space on which these functions are defined as operators so that the commutation 
relations are satisfied. 

In principle, we will get a separate Hilbert space around each geometry from the moduli 
space. Low-lying states in this Hilbert space will not allow any semiclassical interpretation, 
they will be similar to states of a few field quanta (e.g. gravitons) propagating in Minkowski 
space. However, in the limit of large occupation numbers we can construct coherent states, 
which can be interpreted as describing a neighboring classical geometry. In this sense, 
neighboring geometries are contained in each other's Hilbert spaces. This shows that all 
these Hilbert spaces will be isomorphic. 

Let us now discuss the choice of a Hamiltonian operator on the Hilbert space. First 
of all, when we are dealing with a theory of gravity, it is by no means guaranteed that 
a preferred Hamiltonian will exist. For instance, this seems to be the situation for the 
plane gravitational wave case analyzed in |2(:]. For the concept of energy to make sense, all 



spacetimes in the considered moduli space must have a common asymptotic infinity with 
a timelike Killing vector. In this case the Hamiltonian can be defined by computing the 
classical energy of the considered solutions. For the LLM case, such a computation has 
been done in 

When a Hamiltonian is available, the process of quantization can be taken further 
by discussing the energy eigenstates. To avoid possible confusion we would like to stress 
again: the Hamiltonian is an independent piece of information which supplements the 
symplectic form computed by the on-shell quantization method. There are cases when 
it can be defined in the classical theory (and carried over to the quantum theory by the 
correspondence principle), but a separate computation is necessary in order to do that. 

2.10 Corrections 

In which sense does on-shell quantization approximate the true picture attainable in the 
complete theory? And what are the corrections which should be included if one is to go 
beyond this approximation? These are interesting questions, which most likely have to be 
analyzed on a case- by-case basis. Here we will limit ourselves to a few general remarks^''. 

It will be convenient to phrase the discussion in the Hamiltonian language, which may 
not be the most general situation (see Section |2.9D , but is sufficient for the applications we 
have in mind. The full Hamiltonian can be schematically represented as 

F = Fo((/<||) + Fo(0^)+Fint. (2.27) 

Here ifo('/'||) is the quadratic Hamiltonian for the fluctuation modes along the moduli 
space Ai (around a given background); HQ{(j)±) is the quadratic Hamiltonian for the modes 
corresponding to fluctuations orthogonal to Ai; H^^t is the interaction Hamiltonian. 

On the quadratic level there is a complete decoupling between </>|| and (j)^. It is only 
the modes that we are quantizing using the on-shell quantization method, while (f)± are 
^°In the context of minisuperspace approximation to pure gravity, some discussion may also be found in 
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effectively frozen. The result of this quantization will be a Hilbert space 7i\\ and a spectrum 
of energy eigenstates^^. 

It is more or less clear how this picture will have to be modified, if the corrections 
are to be considered. First of all, we will have to introduce a Hilbert space 7i± for the 
orthogonal modes, with its own free spectrum. The full Hilbert space is then the direct 
product TLw TL^. On-shell quantization can be thought of as turning off -fTint and putting 
all (j)x_ in their vacuum state, so that the full state is a product 

|</>l|)® |0x). (2.28) 

Taking i^int into account allows transitions, exciting the orthogonal modes and inducing 
mixings between various states. To be able to treat i^int as a perturbation, a small expan- 
sion parameter should be available. In this case the new energy eigenstates will be small 



perturbations of the original product states ( 2.2S| ), with slightly shifted energy levels. In 



theories of gravity, the role of such an expansion parameter can be played by ^Apiancki the 
characteristic curvature radius d. of the background spacetime measured in Planck units. In 
the string theory context, this will correspond to the a'-expansion of the effective action. 

In supersymmetric situations, the spectrum may be protected and energy shifts should 
not occur. For instance, we know from AdS/CFT that this should be the case for the LLM 
solutions. Indeed, the N = A SYM states dual to these geometries are 1/2 BPS, and 
their energies cannot depend on a continuous parameter. Notice that it would be wrong to 
conclude that the LLM geometries do not get modified once the a'-corrections are taken 
into account — they will, except in the fully supersymmetric cases of AdS^ x and the 
plane wave. It is only the energy eigenstates which should remain protected. It would be 
interesting to show this directly from the gravity side. 

The discussion of corrections becomes increasingly subtle when loop effects are taken 
into account^^. In the AdS/CFT context, the size of these effects is controlled by 1/A^. In 
this paper we mostly discuss the N = oo limit, corresponding to the classical SUGRA. In 
principle, it should be possible to consider ^> 1 case as a perturbation over = oo. This 
would require inclusion of massive string states needed for the UV completion of the theory. 
In the LLM case, we know that this should lead to a reduction in the number of states with 
energies above A^ (see Section \i.3\ ). It would be extremely interesting to understand how 
such a reduction can be achieved in a perturbative treatment, and to reproduce it from the 
gravity side. 



2.11 Summary of on-shell quantization method 

To summarize, we have the following recipe for quantizing a moduli space of solutions 
of any supergravity theory. First of all, we have to find a general expression for the 
symplectic current of the theory. This is computed from the supergravity Lagrangian by 

We do not discuss here the other aspect of the on-shell quantization, namely the correspondence between 
coherent states and classical geometries mentioned in Section |2.9| . 

^^Note that any discussion of such effects is ambiguous and ill-defined in pure gravity, due to its non- 
renormalizability. 
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the general formula ( 2.17] ) and will contain the Crnkovic-Witten gravitational symplectic 



current ( p. 15 ), as well as additional terms for the other fields of the theory. 

Second, we have to evaluate the symplectic current on the moduli space of solutions. 
This is done by expressing the variations of all the fields in the symplectic current via the 
variations of the arbitrary functions describing the moduli space. 

Finally, we have to integrate the symplectic current on a Cauchy surface S to obtain 
the symplectic form, which will be a closed 2-form on the moduli space. This symplectic 
form can then be quantized in the standard way. 

3. Quantization of "Bubbling AdS" 
3.1 Supergravity solutions 

Having set up the general framework in the previous section, we will now apply the on- 
shell quantization method to quantize the "Bubbling AdS" family of supergravity solutions 
found by LLM g. This family includes ah regular 1/2 BPS solutions of Type JIB SUGRA 
with 50(4) X SO {4) x M symmetry. These solutions have constant dilaton and axion and 
vanishing 3-form. The metric has the form: 

ds^ = -h-^{dt + Vidx'f + h'^idy'^ + dx'dx') + ye^dQl + ye'^dQl , (3.1) 

where i = 1,2 , dilg and dCl'^ are the metrics on two unit 3-spheres S^,S^. The functions 
h and G are determined in terms of a single function z{xi,X2,y), y > 0: 

The z and Vi are in turn fixed in terms of one function Z(xi,X2), which can only take the 
values ±1, and is also the boundary value of z on the y = plane. Namely, we have^^: 

vr (x^ + y'^Y ' 
yi = — t 2A2 *^- (3-3) 
Apart from the metric, only the 5-form is turned on: 

= F ^ dVL'i + F ^ dVL?, , 
F = dB, F = dB, (3.4) 

dQ^, dQ-^ being the volume forms of the spheres. The one-forms B, B are defined up to a 
gauge transformation. A convenient choice is the axial gauge By = By = 0, the remaining 
components being then given by 

Bt = --/e^'', Bt = --/e-'"". (3.5) 

„ y'^Vt U, XI ~ y^Vj Ui XI 

4 (i - z) 4 4 ' 4 (1 + 2;) 4 4 ' 



^^We use here the standard notation for two-dimensional convolution: / * g{x) — f f{x — x')g{x')d^x' 
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(see Q), where 



U,^'-^^*Z. (3.7) 



The following linear relations, evident from ( |3.3D and (3/7), will play an important role 
below: 

diz = -ye'^dyVj , 

1 1 

d\iVj] = ^(^iVj - djVi) = ^^ijdyz , (3.8) 

dyUi = -2yVi . 

3.2 Moduli space 

The moduli space is parametrized by collections of droplets of arbitrary shape in the xi,X2 
plane, defined by the condition that Z = —1/2 on the union of all droplets (which we denote 
P) and Z = 1/2 in the remaining part of the plane, 'D'^. We will assume that the boundaries 
of all droplets in D are smooth curves, so that the corresponding solution has a smooth 
geometry We would like to discuss the quantization of this moduli space. As discussed 
in Section there are two aspects to the moduli space quantization. First we must detect 
all nontrivial fluxes. These can be quantized using an appropriate Dirac quantization 
condition. On the other hand, deformations within the moduli space corresponding to 
keeping all the fluxes fixed should be quantized using the on-shell quantization method. 

The topology of the "Bubbling AdS" spacetimes and the corresponding fluxes have 
been already determined in Q. It turns out that for every black or white region in the 
droplet plane there is a non-contractible 5-cycle supporting an flux proportional to the 
area of the corresponding region. For example, for each droplet (i.e. a black region) the 
corresponding 5-manifold is constructed by fibering the sphere over a surface S2 capping 
the droplet as shown in Fig. ^. This 5-manifold is nonsingular, since the shrinks to zero 
size on the y = plane outside of the droplets in an appropriate way. Analogously, the 
5-manifolds corresponding to the white regions are constructed by fibering the over 
surfaces capping these white 'holes' inside the droplets. Using flux quantization, one shows 
that the area of each black or white region must be quantized Q^^: 

Ai = ^^rii, Ui G N. (3.9) 

These integers can be thought of as the numbers of giant gravitons wrapping the 5*^ and 
the S^. The sum of rij's corresponding to the black regions should be equal N, the total 
number of D3 branes making up the configuration. 

We see that there exists an infinite-dimensional class of deformations within the moduli 
space keeping all the fluxes flxed: these are precisely the deformations of the boundaries 
of the droplets which keep fixed the areas of all black and white regions. It is these 
deformations that we will quantize using the on-shell quantization method. 



The relation kio = 8n''/'^£% is useful in comparing some of our equations to 



-15- 



Figure 2: The surface E2 caps the droplet in the xi,X2 plane by extending into > 0. Fibering 
the over S, we get a closed 5- manifold supporting an F5 flux proportional to the area of the 
droplet H). 



3.3 Dual description in terms of free fermions 

The "Bubbling AdS" solutions corresponding to a collection of finite-size droplets are 
asymptotically AdS^ x S^. The AdS/CFT correspondence can be used to relate these 
Type IIB supergravity solutions to states of the J\f = 4 supersymmetric Yang-Mills theory 
on X M. Since the geometries are 1/2 BPS, the dual operators on the Yang-Mills side 
are chiral primaries with conformal weight equal to their U{1) R-charge: A = J. As is 
well known ^ , this sector of = 4 super Yang- Mills admits a very simple description 
as a system of non-relativistic free fermions moving in a harmonic oscillator potential. 
In the large N limit (which is what we are mostly concerned with in this paper) the states 
of the many-fermion system are well described as droplets in the one-particle phase space. 
In fact, these droplets are the same droplets which characterize the gravity solutions. For 
example, it was checked in that the excitation energy of the "Bubbling AdS" solutions 
over AdS^ x is equal to 



A = J ' 



(3.10) 



which is precisely the excitation energy of the fermionic state described by the droplet T) 
over the ground state described by the circular droplet of the same area centered at the 
origin. The one-dimensional Planck constant h here can be fixed by comparing the area 
quantization condition ( p.9[ ) valid on the supergravity side with the semiclassical phase 
space quantization condition A = 27rhn, which should be true on the fermion side. This 
gives 1]: 

In this paper we would like to compare the symplectic structures on the droplet space 
arising from the supergravity and the fermion side. On the supergravity side we will use the 
on-shell quantization method described in Section On the fermion side the corresponding 
symplectic form can be found by using the so-called hydrodynamic approach, commonly 
used to describe edge states in Quantum Hall systems (see e.g. ||3^). The idea of this 
method is to identify the one-fermion description of the classical dynamics of the droplets 
with the collective one. In the one-fermion picture, the droplet motion is described by 
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solving the individual Hamilton equations^^ 

2 2 

p=-dgH, q = dpH, H = ^^. (3.12) 

for the fermions localized close to the droplet boundary 9P. It is convenient to first look 
at the subclass of the droplets whose boundary curve can be described in polar coordinates 
as r = r{(j)) for a single-valued function r((p). For these droplets ( 3.12 ) implies that r((/>) 
evolves in time according to the classical chiral boson equation 

r = h-^r'. (3.13) 

(A cautionary remark: as we will see below, it would be wrong to conclude from this 
that the in quantum theory r(</>) satisfies the chiral boson commutation relations.) In 
the collective pictures, one would like to recover the same equation (|3.13| ) as a Hamilton 
equation of the form 

r = {r,Htot} (3.14) 

where fftot is the total energy of the droplet state, given by the integral of the one-particle 
Hamiltonian: 



Hi 



tot 



dp dq + 



while the Poisson bracket {,} is to be determined from consistency of the two pictures. 
As we showed in (see also |Q, and especially [^, Eq. (4.3)), the appropriate Poisson 
bracket which generates the correct equation is 

{r'^{^),r'^{^)] = Mh5'{<t)-^). (3.16) 

The symplectic form corresponding to these Poisson brackets can be written as^^: 

uJiorm = ^:^'j> j>d^d^S\g■a{ct>-$)5[r\<|))\^5[r^{^ . (3.17) 

The total area of the droplet is proportional to the total number of fermions and must be 
kept fixed. This means that the variations in ( 3.17| ) must satisfy the constraint 

# 5 [r^ ((/>)]= 0. (3.18) 

Under this constraint, the symplectic form ( 3.17| ) is well defined as written, in spite of 
the fact that the kernel Sign(0 — 0) is not periodic. In fact we can add to the kernel any 



^^This normalization of the Hamiltonian has to be chosen so that the energy levels become integer-spaced 
and can be in one-to-one correspondence with the integer dimensions of field theory chiral primaries. 

^^The Poisson brackets are encoded by the symplectic form in the following schematic way: {qi, qj}p.B. ~ 
Aij corresponds to = ^A~^dqi A dqj. To apply this standard rule in our situation, note that the inverse 
of the kernel 5'{(j) — (j>) is ^Sign((^ — 0). 
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function of the form gi{4')+g2i<P) without changing the answer. If desired, one can consider 
adding (cj) — 4>)/Tr, which will make the kernel explicitly periodic. 

It is interesting to note that although the symplectic form ( 3.17| ) is easiest to derive 
for the particular case of the harmonic oscillator one-fermion Hamiltonian, it is in fact 
completely general and will describe the motion of droplets of noninteracting fermions 
described by an arbitrary one-particle Hamiltonian H{p,q) |36|,[||],[37|. 

Let us rewrite (3.17) in a slightly different form, introducing instead of a natural 



parameter s measuring arc length along dT>. From elementary geometry we have 

ds 5r 



(3.19) 



where by 6^±_ we denote the variation of &D in the outside normal direction. Thus ( 3.17] ) 
can be equivalently written as 



t^ferm = ^ j> j> ds ds Sign{s - s) (57_l(s) a 6j±{s) 
The constraint ( 3.1^ ) takes the form 



ds 57j 



(3.20) 



(3.21) 



The advantage of expression ( p. 20 ) is that it is completely general — it makes no reference 
to polar coordinates and can be also used for droplets with multiple-valued r{(j)). It is this 
expression that we will have to reproduce from the supergravity side. 
The Poisson bracket following from ( |3.20| ), 



(3.22) 



is equivalent to the bracket derived by Dhar [37|, Eq. (4.2), using a reparametrization 



invariant description of fermion droplet boundaries. In Appendix A.l we show explicitly 
that ( 3.22| ) generates correct equations of motion of droplet boundary in the general case. 

A few words should be said about the symplectic structure in the situation when several 
droplets are present, or when a single droplet has several boundary components. Assume 
that has B connected components described by closed curves 7(''''(,s). In this case one 
has to introduce a separate field S^^\s) for each boundary. The total symplectic form is 



given by the sum of the forms ( 3.20 ) computed for each boundary component: 



1 ^ 



Svr/i 



(6) 



ds ds Sign(s - s) 57^^ (s) A 57^^ (s) . 



Each Sj^f^ (s) has to satisfy its own constraint 



(b) 



6jf\s)ds = 0. 



(3.23) 



(3.24) 



This means that, semiclassically, different droplet boundaries are completely decoupled, 
and it is impossible for a fermion to move from one boundary to another, even though such 
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a transition would keep the total area of the droplet fixed. Technically, condition ( 3.24 ) 
specifies symplectic sheets in the moduli space of fermion droplets. 

The symplectic form ( p. 23 ) on variations satisfying ( |3.24 ) is all we need for future 
comparison with the supergravity side, since as we discussed in the previous subsection, 
only variations keeping both black and white areas fixed can be quantized by the on-shell 
quantization method, due to the fixed flux restriction. 

Knowing the Poisson brackets, it is trivial to quantize the system. We will do this 
around the circular droplet, the general case being identical. We start by expanding r^(0) 
in Fourier series: 

r2(</.) = J]a„e^"^ = (3.25) 

The zero mode is fixed in terms of the droplet area: 

ao = -. (3.26) 

TT 

The n 7^ modes correspond to the area-preserving deformations. Substituting (|3.2E| ) in 
(3.17), we get Poisson brackets in terms of the corresponding Fourier coefficients: 

{a„, am} = -4:hin6n+m , (3.27) 

which upon quantization become commutators: 

[an, am] = Ah^n5n+m , (3.28) 

The Hilbert space of the theory can be constructed as a bosonic Fock space with a^, n > 
{n < 0) as annihilation (creation) operators. The Hamiltonian can be expressed in terms 
of these operators as 

H = ^ a^nCtn + Const. (3.29) 

n>0 

It is interesting to compare the partition function of this theory 

oo ^ 

^(/^) = n • (3.30) 

n=l 

with the exact partition function for a system of iV < oo fermions in the harmonic oscillator 
potential: 

^ 1 

n=l 

Writing the partition functions in the form 

Z{P) = Y,g{E,)e-P^^ , (3.32) 

we can read off the degeneracy g{Ei) of the energy state Ei. We see that the infinite 
N computation correctly predicts the finite energy levels Ei = 1,2,3,... However, it 
reproduces their degeneracies only up to energies Ei < N, overestimating the degeneracy 
of higher energy levels^^. 

^'^This was first pointed out to us by S. Minwalla. 
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3.4 Symplectic form of Type IIB SUGRA 

The LLM solutions satisfy the equations of motion that can be derived from the following 
consistent truncation of the Type IIB SUGRA action^^: 

•5 = ?V / d'^xV^ (r - 4Fi,^..,^iFI*-*^I) , (3.33) 

We can proceed using this action if we impose the selfduality constraint = on the 
solutions. The presence of this constraint does not modify the underlying symplectic form 
of the theory, which can be computed from the action ( 3.33| ). 



Using the general formulas from Section ^, the symplectic form will be equal to^^ 

^ / dJ^iifc + Jf), (3-34) 

Kin J 



2k 



where Jq and are symplectic currents constructed from the gravity and 5-form parts 
of the Langrangian using ( |2.17 ). In fact Jq is nothing but the Crnkovic-Witten current 



( 2.15 ). To find the 5-form current, we take the potentials Aj^^ (F5 = dA) as our basic 



fields. Applying ( |2.17| ), we get immediately 

4 = -85(^F'I^'-'=*I) A5A|fc,...,,,. (3.35) 

Now that we computed the symplectic current, we can simplify it using the selfduality 
constraint, which can be written as 

where e " is the flat 10-dimensional epsilon symbol. Thus we have 

ji^ ^ -8e'l^'-'=*ll"-"^^lji^^,...^,l A5A|,,...fc,,. (3.37) 
3.5 Regularity condition. Boundary term 

To evaluate the symplectic form, we will use T. = {t = const} in (|1|). This is a natural 
choice, in view of the time-independence of the solutions. With this choice, J* is the only 
needed component of the symplectic current. Because of the 50(4) x SO (4) symmetry, 
integration over the 3-spheres is trivial, and we will have 



uj = - — J- / dy d^x{J^G + Jp) ■ (3-38) 
2^10 J 



There is however a subtlety to be discussed, before we start evaluating this integral. 
Namely, we have to make sure that the field perturbations we are using are regular. More 
precisely, this condition means that 5gmn and JA4 should be regular in a local coordinate 
system in which the background metric gmn is regular (such a coordinate system exists 



^*The notation means that the indices have to be ordered: ii < ... < i„. Thus we have 

F = _F|ij...ij.|da;'^ A ... A da;*^ — ^Fi-^,,,i^dx^'^ A ... A dx'^ . The same ordering is assumed in the summation. 
^^The symplectic form and currents appearing in this section were aheady derived in H. 
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because all LLM geometries are regular). If this condition is violated, we may get a wrong 
result for the symplectic form. In fact, as we will explain below, this was precisely the 
origin of the factor 1/2 discrepancy encountered in |^]. The source of trouble is (5^4, which 
has the form 

5Ai = 5B Adn + 5B Adn (3.39) 

and does not satisfy the regularity condition on the y = plane if the axial gauge ex- 
pressions ( |3.6D are used. In fact, 6B and 6B obtained by perturbing ( |3.6D tend to a finite 
nonzero limit as y ^ 0. On the other hand, the spheres S^{S^), whose volume forms 
appear in (|3.39D , shrink to zero size inside (outside) the droplets. This means that the 
four-form (|3.39| ) is singular as y — > 0. 

To get the right answer for the symplectic form, we must use in ( 3.37] ) a regular gauge 
field perturbation 6A]^^. The existence of such a regular perturbation and its properties 
are discussed in detail in Appendix A. 2 . The result of this analysis is that 5A^^^ can be 



chosen to have the form 

6Al^^ = 5B'''^ Adn + 6B'^^ A dn , (3.40) 

where 5B^^^,6B^°^ are regular functions at y > 0, which have zero y = limits on 'D(T>'^) 
as a consequence of the regularity of 6A]^^: 

5B'''%=o = on P, 

5B"'%=o = on 2?^ (3.41) 

(the precise regularity condition is somewhat stronger, see appendix A.2). 



Since 5A]^^ must be related to dA/^ by a gauge transformation, we must have 

dB'''^ = 6B-dX, dB'""^ = 5B-d\ (3.42) 

for some regular functions A, A. 

Now we are ready to find the 5-form part of the symplectic form. First of all, substi- 
tuting (l3lc| ) into (|3!37|) , we get 



4 = 4e'^''^( 5Bl^^ A 6Fb^ - 6Bf^ A 5^^^ ) . (3.43) 
Using (|3l^ ), this can be expressed as 

4 = 4e"^^( 55, A 5Fbc - 5Ba A 6Fbc ) - ie^^'idaX A 5Fbc - 9a A A 5Fbc ) 

= 8e'^{dy6Bi A 6Bj - 6Bi A dySBj) - Ae^^'^dai A A (5^^^ - A A SF^^ ) , (3.44) 

where we simplified the first term using the axial gauge condition By = By = and used 
the closedness of F and F to rewrite the second term as a total derivative. Contribution 
of this term to the integral in (3.38) will be given by a boundary term at y = 0. More 



precisely, we will have 

[ dyd^xJ'p = I dyd^xJp'^^+ [ d^xjf"^ , (3.45) 
J J Jv=0 
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where 



jbd^-y = 8 (A A 5Fi2 - A A 6Fi2 ) . (3.46) 



If we had used the singular perturbation 6A4 in (3.37), we would have missed the boundary 
term in ( p.45| ). This is what caused the factor 1/2 discrepancy encountered in Section 6 of 
0. 



Eq. ( 3.45 ) is the final general expression for the 5-form part of the symplectic form. 
We see that to evaluate it, we only need to know the functions A (A) at y = 0. Their values 
on P (P'^) can be determined from ( ^.41 ). Using ( 3.42| ), we see that we must solve the 
equations 

5Bi = diX on V, 

6Bi = diX on (3.47) 

In the complementary regions A (A) can be defined arbitrarily, since the field strengths 
multiplying these functions in Jp'^'^^ will vanish there as a consequence of (3.47). This 



is of course not surprising, since making A or A nonzero in the complementary regions 
corresponds to a perfectly regular gauge transformation. 

Let us repeat the logic of the preceding discussion. We observed that the axial-gauge 
perturbation 5A4 is singular at y = 0. We must therefore gauge-transform into a regular 
perturbation 5A]^^. Such a regular form must exist, since all the fluxes are kept fixed. In 
principle, it is 5A]^^ which should be used to compute the symplectic form. However, 
in practice it would be convenient to be able to use the explicit axial gauge expressions. 
Eq. (|3l5|) shows that we can keep using (5^4 to compute the symplectic current, but we 
have to pay a small price of adding a term to the symplectic form localized on the lower- 
dimensional submanifold where 5A4 becomes singular, i.e. at y = 0. We will refer to this 
term as the boundary term, although of course it should be remembered that y = is only 
a coordinate boundary and the spacetime is completely regular there. 

It remains to discuss the regularity of Sg-mn- This perturbation is manifestly regular 
everywhere at y = except maybe at dD. Thus we see that the situation is better than 
with 5A4, already because the singularity occurs on a smaller set. In fact, one can show that 
the remaining singularity on dV is harmless, and does not give rise to any additional terms 
in the symplectic form. The argument, which we will not present here in detail, proceeds 
by constructing a linear gauge transformation which completely removes the singularity in 
^9mm and by showing that this gauge transformation has no effect on the symplectic form. 
Thus the gravitational part of the symplectic form can be computed by using the naive 



expressions for 5gmn in the Crnkovic-Witten symplectic current ( 2.15 ). 
3.6 Wavy line approximation 

We will now use our method to compute the symplectic form around the plane wave 
background, described by the droplet filling the lower half-plane: T) = {x2 < 0}. The 
fluctuations around this background correspond to deforming dD from X2 = into X2 = 
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e(xi). The symplectic form will be an antisymmetric bilinear functional of e(x). The most 
general form of such a functional is 

. = //<i.di-A>,iO.(.)A.(i), = i)^ (3.48) 

Since xi — > xi + const is a symmetry of the plane wave background, the kernel K{x, x) 
will be translationally invariant: K(x,x) = K{x — x). This suggests that we should work 
in the Fourier-transformed basis, in which uj will have a diagonal form^'': 

u= f ^K{p)e{p)Ae{-p), K{p) = K{-p), (3.49) 



eip) = J dxa'P'^six). (3.50) 
We will now determine the precise form of K by an explicit computation. 
3.6.1 Gravitational part 

In this subsection we will evaluate the gravitational part of the symplectic form, which 
is given by the integral of the Crnkovic-Witten current ( 2.15| ). It should be stressed that 



splitting the full symplectic form into the 'gravitational part' and 'five-form part' has no 
physical meaning; in particular, these parts will not be separately invariant under gauge 
transformations or changes of S. It is done here for purely computational purposes. 

The background values of the fields characterizing the plane wave metric are found 
from (l3.3|) for Z = iSignx2 and are given by (r = -y/x^ + y^) 

z = ^, Vi = -^, V2 = 0. (3.51) 
2r 2r 



To find the perturbations of these fields, we have to expand ( |3.3D to the first order in 
e(x). As we noted above, it is natural to Fourier-transform w.r.t. xi. In this basis, the 
perturbations have the following form |2|: 

6z{p) = -y'^^^e-^'^Mp) , 
SViip) = -X2^^^e-^\MP) , 

SV2{p) = '^e~^^Mp) , (3.52) 



5z{p) = 5z{p,X2,y) = / (ixi e*^^Mz(xi, X2, y), etc. 



where 



To find perturbations of the metric components one should also vary Eqs. (|3.2| ) and use 
the above 5z. We do not present the results of these elementary computations here. 

According to ( |3.38| ), we have to compute the time-component of the Crnkovic-Witten 
current ( 2.15| ). The computation is straightforward. The current is a bilinear expression 

^"To avoid proliferation of tildes, we denote Fourier-transformed functions by the same symbol with p as 
an argument. 
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in e{pi), i = 1,2. Here we present the result of integrating this expression in xi, which 
produces a momentum-conserving (5-function pi + P2 = 0: 

2rV + 2r\p\-l + (2rV + ^r\p\ + 3) 



dxiJr 



dp 
2^ 



.2 



^e-2Hpl,(p) A e{-p) . 



(3.53) 

This expression stiU has to be integrated in y > and X2- The dependence on p can be 
scaled out of the integrand by making the change of variables r — > The remaining 

integral is easy to evaluate in polar coordinates, and the final result is: 



dxi dx2 dy Jq 



dp i . . . . 

e{p) A e{-p) 

Ztt p 



(3.54) 



This result was also derived in |^ by using the general expressions for the symplectic 
current, which we will present in Section 3.7.1 below. 

3.6.2 Five- form part 

Perturbations of the gauge field potentials can be analogously found from (|3.6| ) 



SBiip) 



1 



-r\p\ 



4 

4 



[x2\p\ + (1 + \p\r)]e{p) 



e ■ "-' \x2\p\ - (1 + \p\r)] e{p) 



(3.55) 



5B2{p) = {i Signp) 6 Blip) , ^Mp) = "(^ Signp) 6Bi{p) . (3.56) 
As we found in the bulk term in the gauge field symplectic form expression ( |3.45| ) 



vanishes on these perturbations: 

j d^x dy J^"^^ = (plane wave) . 



(3.57) 



Thus the only contribution will come from the boundary term in ( 3.46 ). The field strengths 
at y = are easily found to be 



5Fi2{p) = -p^X2e-'"'\P\e{j)) {x2 > 0), 
5Fi2(p) =p22;2e^2|ple(p) (x2 < 0). 



(3.58) 



The field strength variations vanish in the regions complementary to the given, and thus the 
corresponding variations of the potentials are pure gauge, in agreement with our discussion 
in Section 3.5. Solving ( |3.47| ), we find the gauge transformation parameters: 



\{p) 
\{p) 



.J_qX2\v\ 

ip 



e{p) {x2 < 0), 



-X2\p\ 



4p 



e(p) 



ix2 > 0). 



(3.59) 



With this information it is possible to evaluate the boundary term in ( |3.45D : 



ly=0 



d^xjf'^ 



AA(5Fi2-8 



j/=0,X2<0 



y=0,X2>0 



AA,5Fi2 = / ^-e{p)^e{-p) • (3.60) 
27r p 
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3.6.3 Plane wave symplectic form 



Adding the contributions of ( 3.54| ) and ( 3.60| ) in (3.3J), we get the symplectic form for 



perturbations around the plane wave background: 

Since i/p is the Fourier transform of ^Signx, this can be rewritten in the coordinate space 
as 

Lo = I I dx dxSignix — x) eix) f\ eix). (3.62) 
^10 J J 

We would now like to compare this answer with the fermion symplectic form Wfcrm given by 



( 3.17 ). These symplectic forms encode Poisson brackets between different Fourier harmon- 
ics of the boundary curve perturbation 7_l(s) = e{x). Strictly speaking, we must compare 
commutation relations, which differ from the Poisson brackets by an extra Planck constant 
factor, hiQ on the supergravity side and h on the fermion side. Since the 10-dimensional 
Planck constant satisfies the usual convention Hiq = 1, the precise relation which should 
be satisfied is 

UJ = fr^UJ{c^ja ■ (3.63) 
Using the value of h from ( |3.11| ), we see that this is indeed true. 

3.7 General droplet case 

In this section we will compute the symplectic form on the full moduli space. 
3.7.1 Current expressions 

The first step is to evaluate the symplectic current as a function of variations of the var- 
ious fields describing the solutions. The gravitational part of the symplectic current was 
evaluated in 0] and is given by 



Jh= y' 



^diVidih-"^) A 6{h^) + 3S{VidiG) A6G + 25{h-% %Fj]) A 8Vi 



A5{di\nh) A 5Vi 



(3.64) 



The last term can be represented as a total derivative by using diVi = 0. In the remaining 
terms we express h and G via z from (|3.2D . We also choose to use the linear relations (^^) 
to eliminate diZ and in favor of dyVi and dyZ. The resulting simplified expression has 
the form: 



Jg = ^ 7T^^J HVj dyVi) ASz- y% S 1 A 5Vj - 4di {5\nhA SVi) ■ (3.65) 
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To compute the five-form symplectic current, it is convenient to rewrite the gauge potential 
perturbations (3.6) in the form 



SBi = - 
ai = 6 

hi = 5 



y^zVi 



(3.66) 



The bulk integrand in (|3.45|) then takes the form: 



jbulk 



Eijitti A dybj - dytti A bj) 
-2eij dytti A bj + dy {eijUi A bj) 



-2eij 6 



9, 



'HI 



2yV^ 



A6 



+ dy {sijUi Abj) , 



(3.67) 



where we used ( |3.8| ) to express dyUi via Vi in the first term of this expression, so that it 
became a function of Vi and z only. 

3.7.2 Total derivative representation 

To find the symplectic form, we will have to integrate the bulk currents Jq and Jp^^^ found 
in the previous section. We see that these currents are given by complicated nonlinear 
expressions, so that performing this integration seems to be a daunting task. Our only 
hope to do the integral is to look for a total derivative representation of the total bulk 
current: 



J 



bulk 



jt I rbulk 



dyP + diP 



(3.68) 



If such a representation exists, then the bulk integral will localize to the (coordinate) 
boundary and will be computable. 

As we will explain now, there is a systematic way to look for representations of this 
kind. First of all, notice that the last terms of ( p.65| ) and (3.67) already have the required 
total derivative form. Let J^^^^ denote the sum of the remaining terms. We can expand 
jbuik ^j^g basis of 2-forms built out of 5z, 5Vi and their y-derivatives: 



J 



bulk 



AeijbVi A 5Vj + BeijdySVi A 5Vj + M,6V, A dy5z + NidydVi A5z + PidVi A 5z. (3.69) 



The coefficients in this expansion are easy to find from ( p.65| ), ( p.67| ), and we won't present 
them here. 

Now we consider the following ansatz: 



fy = a EijbVi A bVj + i3i 6Vi A 6z , 
r = F6ViA6z. 



(3.70) 
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This ansatz is natural if we assume that I'^ is given by local expressions in field variations. 
Since should contain one derivative less than J^""^, it should be built entirely from 5z 
and 6Vi. Comparing the coefficients, we get the following set of equations: 

A = dytt 

B = 2a-yF 

M, = Ni = Pi 

Pi = dyl3^ + d^F (3.71) 



(to get the second equation, one has to express diZ via dyVi using ( |3.8| )). It turns out that 
these equations indeed admit a solution. Namely, from the first three equations one finds: 



a 



A = -^-^^^^, (3.72) 

4 ^ 

One then can check, using ( |3.8D , that the remaining fourth equation is also satisfied, and 
thus we have constructed a total derivative representation for J^"^'^. The total derivative 
representation for J*^""^ follows when we take the last terms in (3.65) and (3.67) into 
account: 

P = EijSVi A 5{aVj) + Eijtti A bj , (3.73) 
P = F5Vi^5z -A5\nh^ 6Vi . 

Here it was possible to simplify the P by noticing that the found coefficients a and (3i are 
in fact related by 

da 
dz 

We end this subsection with the following remark. One could think that since the 
symplectic current is conserved, it can always be represented as a total derivative, i.e. we 
will always have 

*J = dQ. (3.75) 



A = —e^,V, . (3.74) 



This is of course true, at least locally. But there is a crucial point, which makes looking for 
a representation of this kind rather useless in general. Namely, while J of course is given 
by an expression which is local in the variations of the fields, G in general will be non-local. 
For instance, in the chiral boson example considered in Section |2.2| , the symplectic current 
restricted to the chiral boson moduli space has the form 

r = 6f{u) A 5f{u) . (3.76) 
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It is obvious that this cannot be represented as a u-derivative of an expression local in 5f. 

It is thus clear that something special is going on in our case, since we did find a 
total derivative representation local in the variations of the fields. It would be extremely 
interesting to provide an explanation for this phenomenon, and in particular to find out if 
supersymmetry plays any role in it. 

3.7.3 Symplectic form 

Using the total derivative representation found in the previous section, we find that the 
integral of J^'^^^ localizes to the (coordinate) boundary. Since 6Vi decays at x,y — > oo, 
the only contribution will come from 1^ and it is located at y = 0. Adding the Jp'^'^^ 
contribution (see ( p.46[) ), we get 



I dy d^xiJ'c + 4) = [ d\ i-P + J^/'^) 

J Jy=0 



EijdVi A 5{aVj) - Eijai Abj + 8 (A A dFu - A A SFu ) 



(3.77) 



ly=0 

We now proceed with further simplifications, beginning with the last term. We have 



(A A 6Fi2 - A A 6Fi2 ) 



y=o 



A A e'^ di6Bi 



V 



A A e'^ di6Bi 



DC 



e'^ diX A 6Bj + / e'^ diX A SBj + / (A A SB' + XA6B- 

V Jv JdV 



y=o 



(2z) e^^B^ A 6Bj + / (A A 6B- + A A 6B^ 



(3.78) 



Here in the first line we used the fact, following from the general discussion in Section 3.5 



that 5Fi2{5Fi2) will vanish in 'D(T>'^). In the second line we used Stokes' theorem, and in 



the third line we used ( 3.47 ) and also united the integrals over D and D'^ by inserting the 
2z factor to account for the sign difference. The superscript +(— ) on 5B and 5B indicates 
that the limiting boundary value on dD should be taken from outside (inside) P. 

To simplify further, we will need the following expressions, which are easy to compute 
from ( |3.66| ): 



EijUi A bj = £ij5Vi A 5 



+ EijSUi A 5 



4 ^ 



(3.79) 



16z e'^bBi A bB^ = EijdVi A 5 



+ z£ij5Ui A 5 



4 ^ 



6Vi A 5z + zEijSUi A 6Uj . (3.80) 



Plugging these into (|3.77] ), we get the following simplified expression for loq: 
= 



y'^Eii Vj 



^ 7 {6Ui + y^6Vi) A5z+ zSijSUi A 5Uj + 8 (A A 6B- + A A SB+). 

y=0 1- Z Jy=0 JdV 



(3.81) 
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What happened here is that the first terms of ( ^.79 ), ( ^.801) have completely cancelled 
against 1^. The second terms of ( p. 79 ), ( 3.8C| ) also partially cancelled between themselves. 



The first term in (3.81) still looks rather complicated, however it turns out that it 
vanishes. Roughly, this happens because Sz for y ^ tends to a (5-function supported on 
dT>. At the same time the expression multiplying 6z tends to a function which vanishes on 
dT), and thus the integral of the product tends to zero. A careful demonstration of this 
fact is given in Appendix |A.3 . 



Now that we are left with just two terms, let us take a close look at the gauge param- 
eters A, A. Recall that they have to be determined on V{V'^) by solving ( 3.47| ). It is easy 
to see that the first terms of the axial gauge expressions (^]^) vanish in the corresponding 
domains, and thus we have: 



Thus ( 3.47| ) takes the form 



5Bi 
6Bi 



1 



diX 



on V , 
on v. 

on V , 
on V . 



(3.82) 



(3.83) 



As we explain in Appendix A. 4, 6U\y=o is closed and, moreover, has zero integral on all 
topologically nontrivial cycles in V and V^. Thus equations ( 3.83| ) always have a solution. 
Using (3.83), the second term in (|3.81|) can be transformed into an integral over dV: 



zEijSUi A6Uj = - 



EijSUi A 5Uj - - 



V 



EijSUi A 5Uj 



SijdiX A 6Uj + 2 



SijdiX A 6Uj 



+ XA6U~ 

Thus we are reduced to the following representation: 



2 / {XA6U~ 



dV 



X A {26U- + 86B-) + A A {26U^ + 86B^ 



(3.84) 



(3.85) 



As we explain in Appendix |A.5| , 5B and SB are continuous across dD. Using this informa- 
tion and (|3.82| ), we have 

1 



dB+ 
6B- 



6B- 
6B- 



-6U' 



5U+ 



(3. 



On the other hand, the tangential component 5U^^ has a jump discontinuity on dT>, the 
opposite-side limits being given by (see Appendix I 



(3.87) 
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In view of ( |3.83|) , we can express X,X\dv via SU\\\qx> as follows 



21 



X{s) 



ds6U^{s) 
ds6U^{s) 



ds (^7^(5) 
ds5^±{s) . 



(3.88) 



Eqs. ( 3.86| ), (|3.87] ), ( |3.88| ) contain all information necessary to evaluate ( 3.85| ) in terms 
of (57_L. We obtain: 

uq = 2 j) ds 57_L (s) A ds 5'~i±_ (s)^ = 2 j) j) ds ds 9{s — s) 5^±_ {s) A (57_l (s) 

ds ds Sign(s - s) (57_l (s) A (57_l (s) , (3.89) 



where we took advantage of the constraint ( 3.18| ) , which allows us to add a constant to the 
kernel without changing the symplectic form. 

We see that loq has exactly the same functional dependence as the fermion symplectic 
form Wferm from Section The supergravity symplectic form uj differs from lvq by a 
constant factor (see ( |3.38| ), ( |3.77| )). As discussed in Section 3.6.S, the relative normalization 
of UJ and cjferm should also match. It is easy to see that the precise matching condition as 



expressed by Eq. ( |3.63D is indeed satisfied, just like in the plane wave case. 

To close the discussion, we would like to note that Eq. ( |3.85| ) is completely general and 
in particular applies when the droplet boundary has several disconnected components. In 
this case we can still use Eq. ( p.86 )-( |3^ ) on each boundary component separately, and 
reproduce the general fermionic symplectic form ( p. 23 ). 



4. Conclusions and future directions 



In this paper we have laid out a general procedure for quantizing moduli spaces of solutions 
of supergravity, on which all fluxes are kept fixed. This procedure, the 'on-shell quantiza- 
tion', entails writing down the CWZ symplectic current for the SUGRA action, evaluating 
it on the moduli space and integrating it over a hypersurface to obtain the symplectic form 
on the moduli space. The Poisson brackets which are obtained from this symplectic form 
are then promoted to commutators, which are used to build the whole Hilbert space. 

We have stressed that in order to obtain a non-degenerate symplectic form, it is im- 
portant that the solutions be dynamical (i.e. not static). We also demonstrated that the 
on-shell quantization is really equivalent to the more conventional method of expanding the 
action to second order around a solution, but has the advantage that it automatically gives 
the truncation of all fluctuation modes into the moduli space of interest. We discussed 
the corrections to this quantization, appearing upon introducing the full set of possible 
fluctuations including those transverse to the moduli space, and explained that although 

^^The choice of the lower hmit of these integrals is unimportant and does not affect the value of the 
symplectic form below. 



-30- 



in general energy levels could shift, the spectrum may be protected against such corrections 
in presence of supersymmetry. 

As an application of the developed procedure, we analyzed the family of "Bubbling 
AdS" SUGRA solutions found by LLM Q. We noted that one must work in a gauge choice 
where the variations of the metric and RR gauge fields are regular, and wrote down the 
expressions for the symplectic current taking this into account. Remarkably, it turned out 
the symplectic current around any droplet or collection of droplets can be written as a total 
derivative of a local expression. This enabled us to integrate it and find the symplectic 



form in full generality, Eq. ( 3.89 ). Comparing this to what we expect from the fermion 
picture, we find exact agreement. 

In other words, we have proved that the Hilbert space around any point in the moduli 
space of LLM solutions is the same as that of free fermions in a harmonic oscillator, both 
quantized around the same droplet configuration. We stress again that we only used the 
SUGRA action and solutions, and in no way used AdS/CFT or string theory to obtain this 
result. Thus one can view this as a derivation of AdS/CFT in the N ^ oo limit, in this 
restricted subsector of 1/2 BPS states^^. 

In this sense, our on-shell quantization is a powerful tool, especially in cases where a 
dual description of a SUGRA system is not known. It some situations it can be used to 
derive such a dual picture, where the microscopic degrees of freedom are manifest. It is also 
very useful as it enables one to count the dimensions of certain subsectors of the Hilbert 
space — certain ensembles of microscopic states with common macroscopic characteristics 
— and thus deduce their entropy. This is of course a very interesting question, especially 
when one deals with black hole geometries, such as the 2-charge D1-D5 and 3-charge Dl- 



D5-P black holes that we have discussed in section U^. In future work, we intend to apply 



this method to quantize the 2-charge D1-D5 system [38|. It would be nice to explore other 
systems and setups where the on-shell quantization method could be applied. In particular 
it would be interesting to understand the subtleties in applying our method to spacetimes 
with horizons and to spacetimes which include both electric and magnetic sources'^^. 

Another very interesting question to explore is the 1/A^ corrections to our results. 
As we mentioned in Sections ^.10| and |3.3| , such corrections should effectively reduce the 
dimension of the Hilbert space by decreasing degeneracies of the energy levels above A^. 
It would be extremely interesting to understand how such a reduction may occur in a 
perturbative fashion'^^. 
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A. Appendix 

A.l Fermion symplectic form for multiple- valued r((j)) 

Let us show that the Poisson bracket ( 3. 22] ) generates the correct equations of motion of 



the droplet boundary for multiple- valued r{(j)). The energy of such droplets is still given 
by Eq. (|3.15 ), which should be understood as a line integral along the boundary, so that 



dcj) = {d(j)/ds)ds, and d(j)/ds changes sign at points separating different branches of r{(j)). 
Taking this change of sign into account, the energy can be rewritten in terms of separate 
branches r = ri{(j)) as 

^tot = Y^T.'^ [ ^'^^'('^)' = Sign^ , (A.l) 

where li are the intervals on which the corresponding branches are defined. Furthermore, 



each branch satisfies the same equation of motion ( 3.13 ) as before: 

r- = fr^r'i. (A.2) 

To rewrite ( p.22| ) in terms of r^, note that the relation ( p. 19] ) remains true for each 
branch (when s crosses from one branch to the next, both d(j)/ds and 6r/5^± in ( 3.1S| ) 
change sign). The 5' transforms as 

^'(.(0) - sm = - 4>) , (A.3) 

s (0)s (</>) 

and we get: 

{6rH4>),Sr^{4>)} = 87Th6ijei6'{(t> - 4>) • (A.4) 



The Hamilton equations corresponding to this bracket and (A.l) indeed coincide with 

(Ia: 



A.2 Regular gauge field perturbations 

In this section we demonstrate the existence and properties of a regular gauge potential 
perturbation 6A]^^ used in Section 3^. First, we note that the closed 5-form 6F is regular, 



since it is a difference between two regular 5-forms corresponding to neighboring LLM 
geometries. Moreover, since we assume that all the fluxes are kept fixed, 6F has zero 
integral on all closed cycles. Thus by de Rahm's theorem, it is in the trivial cohomology 
class: there exists a regular 4- form 6A]^^ so that 

6F = d{6Al^^) . (A.5) 

Furthermore, this 4-form may be assumed to be SO{A) x 50(4) invariant, and thus of 
the form ( 3.4C| ). This can be achieved by averaging over this group. More precisely, every 
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group element g G SO (4) x 50(4) acts in an obvious fashion on the LLM manifold, rotating 
the spheres and S^. This action leaves 6F invariant: g*{6F) = 6F. Since the pull-back 
action commutes with exterior differentiation, we get 

5F = d{g*{SAl'<^)) . (A.6) 

Averaging this last equation over g, we get the desired SO (A) x S'0(4)-invariant represen- 
tation. 

Finally, let us discuss the requirements imposed on the components of 5B^^^, dB"^^^ by 
the regularity of ( 3.40 ). Consider the case of 5B'^'^^ first. Although the sphere 5^ shrinks 
to zero size on P, it does so in a regular fashion, uniting with the y coordinate to form a 
metric which looks locally like the origin of : 

c{x){dy^ + y^d^l). (A.7) 

The volume form of this R^, which is of course a regular 4-form, looks in these coordinates 
like y'^dy A dVL^^. Comparing this with ( ^.40 ), we see that in order for 5A^^^ to be regular 
as y — > 0, 6By^^ should decay in this limit at least just as fast: 

To get a condition for the decay of 5Bl^^, dBl*^^, we should find a minimal power of y which, 
multiplying d^l^, turns it into a regular 3- form (notice that the h stays finite as y ^ 0, and 
the metric in the t and directions is completely regular). In terms of the Cartesian 
coordinates Xf^, /x = 1, . . . , 4, we have 

dHg ~ {X^)-^e^''^''X^ dXy A dXx A dX„ . (A.9) 

Thus it is clear that such a minimal regular multiple is y'^dQ^, and the remaining regularity 
condition is: 

SBl'^dBl"^ = 0{y^) on P. 
The case of 5B^^^ is completely analogous, and the regularity conditions are 

,5B[•^^5B^*^ = 0(/), (55^°s = 0(/) onP^ (A.IO) 

In fact, only a weak form ( p.41| ) of the above conditions were used in Section 3.5. 

A. 3 A vanishing integral 

In this section we show that the first term in ( |3.81| ) vanishes: 

ly=0 i 



f y^EilYl (^5u, + y^6Vi) A5z = (A.ll) 

Jy=0 A ~ ^ 



The integral here should be understood as the y — > limit of an integral over the plane at 
y = y. We will assume that P is a collection of finite-size droplets, and that the boundary 
dD consists of one or more smooth closed curves. We will carry out the analysis separately 
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Figure 3: For (a collection of) droplets of characteristic size I, we split the plane into three regions: 
the near-infinity region Ri ~ {x : dist(x,I?) 3> the near-9I? region R2 ~ {x : dist(a;, dT>) <C ^]^ 
and the intermediate region i?3 covering the rest of the plane. To show that integral ( |A.ll ) vanishes, 
we use large- a; asymptotics in wavy line approximation in R2, and uniform convergence in i?3. 



for the following three regions of the xi,X2 plane, depending on the distance from dV: 
near-infinity region Ri, neai-dV region R2, and an intermediate region (see Fig. |3|). 

The region i?3 is the simplest, since here 5z = 0{y^) tends uniformly to zero, while 
the remaining factors tend uniformly to a finite limit. Thus the integral over this bounded 
region vanishes as 0{y^). 

The region Ri consists of points x with much larger than the characteristic size of 
T). In this case a more careful analysis is needed, since we must check that the integral 
converges near infinity. First we estimate z by 

1 1 y2 1 1 y2 

^=0 T^~^ — 2y2*^'^<^~~A 1* XV, A.12 

2 TT (x^ -|- y^j^ 2 47r x^ 

where xv is the characteristic function, equal 1 on I? and outside. This estimate is true 
as long as y is smaller than dist(x, 2?), and thus applies in Ri for y i- From this estimate 
we get 



2 / 1 \ -1 



XV ] =0{x^). (A.13) 



Further, we have 



\n< ^y^\,y *xv=o{\xr'), 

\Ui\<^^*Xv =0{\x\-^). (A.14) 
x^ + 

Finally, from ( |Al^ ) and ( [O^ ) we can also conclude that 

6Vi = 0{\x\-^) , 5Ui = 0{\x\-'^), 5z = 0{y'^\x\-^). (A.15) 

Combining all these estimates, we see that the integrand in ([A.ll| ) is 0{y'^\x\~'^). Thus the 
integral over Ri is convergent and tends to zero as 0{y^). 

Finally consider the region i?2, consisting of points x which are much closer to cfD than 
the droplet size. In this region we can neglect the curvature of the boundary. Without loss 
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of generality, we will orient the axes so that the droplet covers the lower half-plane (up to 
higher-order terms). The situation now becomes identical to the wavy line approximation 
used in Section |3.6| . Borrowing results from that section, we obtain the following estimates 
for various factors in ( [A.ll ): 



V 



4 ^ 



5V, = 0(r-2) , 
5z = 0{y'^r-^) 



V, = 0{r-^) , 
= 0{1) , 



(A.16) 



For example, the last of these estimates is obtained from ( 3.52| ) as follows: 



|(5z(xi,X2, 



y 



dp 
2^ 



e'^P^^{l + r\p\)e-''\P\e{p) 



< max 
2r^ V 



(1 + r\p\) e 



-r\p\ 



dp 
2^ 



(A.17) 



which gives us the required estimate under the natural assumption that the perturbation 
e(x) is smooth, so that the Fourier tranform e{p) is integrable. The estimates for 5Vi and 
6Ui can be proved in exactly the same way. In the case of 5Ui one should use the following 
Fourier transform expressions (see P] , Appendix B) 



5Ui{p) 



Ple(p) 



'5f^2(p) = *(Signp)e '''^'e(p) . 



(A.18) 



Using ([Al6| ), we see that the integrand of ( |A.llD is 0{y'^r ^). This gives the following 
upper bound for the integral over R2: 



ds 



dv 



dx±- 



0{y), 



(A.19) 



and we conclude that this part of ( A.ll| ) also vanishes as y — > 0. The proof of ( |A.ll ) is 
now complete. 



A. 4 Existence of solutions for A and A 

In this section we prove that ( ^.83 ) always admits a solution. To show this, we use the 
following explicit form of 5U\y=Q, which follows directly from 



B 



IT ^— • 

Ki{x) = 



EijXj 



Ki[x -J 



ib)i 



Sl1\s)ds.^ 



(A.20) 



Here we use the same notation as in ( ^.23 ) to parametrize dT? and its pertubations, which 
are as usual assumed to satisfy the constraints ( 3.24| ). 
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Figure 4: The vector field SU has vanishing integrals over topologically trivial (ri^2) as well as 
nontrivial (Fa) contours in the complement of dV. This is enough to show that ( |3.83| ) has a global 
solution. 



Now we observe that the convolution kernel in ( |A.2C| ) is closed away from the origin: 



(A.21) 



From this it follows that the integral of K over a closed contour T depends only on its 
winding number n around the origin: 



Using this in ( A.2[l| ), we get 



K = 2Tm . 



6U = 



(A.22) 



(A.23) 



for any closed contour lying entirely in D or P^. Indeed, such a contour will have definite 
winding numbers n^^^ around every component of dT?. The integral in ( |A.23| ) will be given 



by a linear combination of constraints (3.24) with these winding numbers as coefficients, 
and thus vanishes. 

Eq. ( A.23| ) applied to arbitrary topologically trivial contours F implies that 5U is closed, 
which guarantees that ( 3.83| ) has a solution locally. Applied to topologically nontrivial 
contours, ( A.23| ) shows that this local solution has trivial monodromy and is thus global, 
which concludes the argument (see Fig. |^. 

A. 5 Limits on dV 



It is obvious from (|3.6D , (|3.7|) that B,B, U stay finite when y — > 0, at least when x ^ dD. 
Consequently, their perturbations 5B,6B,5U will also have finite y = limits. In this 
section we study how these latter limiting values behave on dD. Similarly to how we 



analyzed region i?3 in Appendix A. 3, we can neglect the curvature of the boundary in the 
limit X &D, and use the wavy line approximation of Section |3.6| . 

The Fourier transforms of 5B.i and 5Bi with respect to the longitudinal variable xi 
are given by ( |3.56| ), where we have to put y = 0. These Fourier transforms are manifestly 
continuous at X2 = 0. Thus we conclude that 5Bi and 5Bi are continuous across dV. 

The Fourier transforms of 6Ui are given by ( A. 18 ). Setting y = 0, we see that 5U2 is 
continuous at X2 = 0, while 5Ui has a jump discontinuity: 

6U^ = ±e. (A. 24) 
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This is in agreement with ( ^.88 ), since SU\\ = —SUi corresponds to the positive orientation 
of the droplet boundary. 
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